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Abstract 

The stress-strain relationship of biological soft tissues affected by 
Marfan's syndrome is believed to be non-convex. More specifically, 
Haughton and Merodio recently proposed a strain-energy density lead- 
ing to localized strain softening, in order to model the unusual me- 
chanical behavior of these isotropic, incompressible tissues. Here we 
investigate how this choice of strain energy affects the results of some 
instabilities studies, such as those concerned with the compression of 
infinite and semi-infinite solids, slabs, and cylinders, or with the bend- 
ing of blocks, and draw comparisons with known results established 
previously for the case of a classical neo-Hookean solid. We find that 
the localized strain softening effect leads to early instability only when 
instability occurs at severe compression ratios for neo-Hookean solids, 
as is the case for bulk, surface, and bending instabilities. 
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1 Introduction 



The Marfan pathology is a connective tissue disorder caused by a defect in 
the fibrillin glicoprotein, a major and essential component of all connective 
tissues. As a result, the arterial waLSS of humans with Marfan syndrome 
exhibit weaker mechanical properties and behavior compared to healthy ar- 
terial waLSS. Dilatation of the arteries may further weaken the inner part of 
the artery wall and may eventually cause tearing and even sudden rupture of 
the artery. It is well documented that because of connective tissue changes, 
individuals with Marfan syndrome have high risk of aortic dilatation and 
dissection caused by increased blood pressure and heart rate [Mizuguchi and 
Matsumoto, 2007]. Other problems occurring in the aorta of patients with 
Marfan syndrome include the development of aneurysms, where the artery 
wall thickens and the artery diameter increases [Watton et al, 2004; Watton 
and Hill, 2009] . 

Preliminary results on biological soft tissues affected by Marfan syndrome 
reveal striking differences with healthy tissues, including (i) an isotropic me- 
chanical behavior and (ii) non-monotone uniaxial stress-strain curves with 
a double curvature. This latter phenomenon can be interpreted as a limited 
amount of strain-stiffening followed by strain-softening. It has been called 
a localized strain softening" behavior, denoted with the abbreviation LSS. In 
contrast, healthy arterial wall tissues are typically modeled with just strain- 
stiffening behavior [Holzapfel et al., 2000] and have a highly anisotropic char- 
acter, due to fiber "reinforcement" along the directions of oriented collagen 
fiber bundles. 

Haughton and Merodio [2009] and Merodio and Haughton [2010] recently 
proposed a strain-energy density leading to localized strain softening, in order 
to model the unusual mechanical behavior of tissues with Marfan syndrome. 
The influence of localized strain softening on the bifurcation of inflated cylin- 
ders under axial loading was studied and related to aneurysm formation, with 
particular reference to the mechanical response of arteries weakened by Mar- 
fan syndrome. It was shown that LSS materials trigger bulging solutions. 
Haughton and Merodio's strain-energy density reads 

where fj, is the initial shear modulus, v > is the localized strain softening 
(LSS) parameter, and Ai, A2, A3 are the principal stretch ratios (the square 
roots of the eigenvalues of the Cauchy-Green strain tensors). Note that 
this strain energy density reduces to the neo-Hookean strain energy, W = 
+ A2 + A§ — 3)/2, when v — 0. We also remark that some recent 
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experimental stress-strain tensile curves obtained for mouse models of Marfan 
syndrome do not seem to exhibit LSS [Eberth et ai, 2009]. In effect, the 
precise mechanical properties of tissues affected by Marfan's syndrome are 
still a current topic of intense investigation. Here we simply focus on the 
stability analysis of materials which may be modeled by the strain energy 

0- 

To illustrate how LSS affects the mechanical behavior of a solid with 
strain-energy (jTl), compared to that of a neo-Hookean solid, we consider 
some common types of large static deformations. We call (xi, x 2 , ^3) the 
rectangular coordinate system aligned with the principal axes of the Cauchy- 
Green strain ellipsoid, and we examine in turn how the solid reacts to a plane 
strain deformation in the Xi-direction, with principal stretches 

Ai = A, A 2 = A"\ A 3 = l, (2) 

say, then to an equi-biaxial deformation in the x±- and 23-directions, with 
principal stretches 

Ai = A, A 2 = A" 2 , A 3 = A, (3) 

and finally, to a uni-axial deformation in the Xi-direction, with principal 
stretches 

A 1 = A, A 2 = A^ 2 , A 3 = A- 1 / 2 . (4) 

Figure [I] shows the stress response in the cases g and Q, with plots of 
the principal Cauchy stress a = o~\ as a function of the stretch A = Ai. 
Clearly, the LSS parameter is an influential factor of loss of convexity for 
these deformations, both in compression and in tension, as it leads to the 
apparition of small windows of strain-softening. 

Our objective here is to study the stability of solids with strain energy Q 
when they are subject to large compressions, by considering the possibility 
of superposed small-amplitude (i.e. incremental) static solutions. This is the 
so-called Euler criterion of stability [Beatty, 1996] . We restrict ourselves to 
two-dimensional incremental deformations, for which the associated mechan- 
ical displacement u depends on two spatial coordinates only. With respect 
to stability analysis, three quantities a, (3, 7 play an important role. They 
are given in general by [Dowaikh and Ogden, 1990] 

a = X 2 1 (X 1 W 1 -X 2 W 2 )/(X 2 1 -X 2 2 ), 7 = aAr 2 A 2 , (5) 
2/3 = (WnXj + X 2 2 W 22 ) /2 - X^W^ + XtWt - a, 

where W { = dW/dXi, = d 2 W/dX,iXj = 1,2). These quantities 
have the dimensions of elastic stiffnesses and play a role similar to that of 
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Figure 1: Effect of localized strain softening for the large deformations of 
solids with strain-energy density given by ((!]). Stress-stretch curves in com- 
pression (A < 1) and in tension (A > 1), for plane strain deformation (figure 
on the left) and uni-axial deformation (figure on the right), with v — 
(neo-Hookean solid), and v = 0.25, 0.75, 1.25, 1.75. 



elastic constants in an anisotropic material with orthorhombic symmetry. 
The difference here is that the anisotropy is strain-induced, not intrinsic and 
fixed, so that the values of a, (3, 7 change as the stretches Ai, A2, A3 change 
under increasing compression. For the strain energy 0, we find that they 
are given by 



a 
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fiXl 



1 + 2v{h - 3)e-" (/l - 3)2 
1 + 2v{h - 3)e- /{h - 3)2 



(6) 



2/3 = a + 7 + Aiiv{X\ - X\f [l - 2v(I x - 3) 2 ] e 



-Hh-3) 2 



with I 1 = X\ + X 2 2 + A|. 

In Section [2j we study the bulk stability of an infinite solid with strain- 
energy ([T]) under compression. We find that, in contrast to the neo-Hookean 
solid, such a solid can become internally unstable at sufficiently high com- 
pression ratios, at least as long as v is away from 0.0 (roughly, for v > 0.5). 
Then, in Section |3j we turn to the simplest boundary value problem, with the 
stability analysis of a semi-infinite solid under compression. As expected, we 
find that surface instability occurs earlier than bulk stability, at critical com- 
pression ratios which are above those of a neo-Hookean solid (again, when 
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v is away from 0, say v > 0.1). These two sections suggest that solids with 
strain energy (flj are much more unstable in compression than classical hy- 
perelastic solids. However, the picture is more nuanced when we turn our 
attention to problems of compressive instability in solids with finite size. 

For the compression of rectangular slabs (Section ^ or of cylindrical cir- 
cular tubes (Section [5]), we find hardly any difference with the neo-Hookean 
case, at least for slender bodies with dimensions compatible with the geom- 
etry of biological tissues. Only for the bending instability problem (Section 
[6]) do we find a significant difference with the neo-Hookean case. This can be 
explained by the known close relation between bending and surface stability 
[Gent and Cho, 1999], thus making the investigation of Section [3] worthwhile 
a posteriori. For those finite size stability prototypes, we use dimensions 
which are representative of human arterial tissue dimensions. 

We focus on compression stability issues because arteries are often com- 
pressed, for instance when they contract or when they are bent at tortuous 
sites. In fact, the compression of the right coronary artery by a thoracic 
aneurysm of the ascending aorta has been reported for a patient with Mar- 
fan syndrome [Minami et al, 2007]. The separate problem of stability under 
large (lumen) pressure is treated elsewhere, see for example Merodio and 
Haughton [2010] or Han [2007]. 

The tentative conclusion of this paper is that for an instability occurring 
at low levels of compression in a neo-Hookean solid (such as for the com- 
pression of slender bodies, Sections [4] and [5]) , the LSS effect does not play 
an influential role. In contrast, it makes a solid with strain-energy density 
([I]) more unstable than neo-Hookean solids for instabilities occurring at high 
levels of compression (Sections 



2 Bulk stability 

To investigate bulk stability (also called internal stability by Biot [1965]), we 
study the ellipticity, or lack thereof, of the differential equations governing the 
equilibrium of small-amplitude two-dimensional static solutions, existing in 
the neighborhood of a large compressive strain of a solid of infinite extend. 
Dowaikh and Ogden [1990] show that those equations are strongly elliptic 
when 

a > 0, /3 + v /«7>0. (7) 

Clearly, in view of the expression (J6]) 1 for a, the first inequality is always 
satisfied for solids with strain energy ([I]), as is the case for neo-Hookean 
solids. However, in contrast to those latter solids, the other inequality may 
or may not be satisfied, depending on the magnitude of the LSS parameter 
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Figure 2: Loss of ellipticity (internal instability) for a compressed infinite 
solid with strain-energy density given by Q. For a given stretch ratio A of 
uni-axial compression, plane strain compression, or equi-biaxial compression, 
and for a solid characterized by a given u, the incremental equations of 
elastostatics loose their ellipticity when a point with coordinates (A, v) is 
situated below the thick curves. For the neo-Hookean solid [y = 0), there is 
no loss of ellipticity. 




v. To illustrate this possibility, we specialize the analysis to certain common 
large deformations. 

First, take the large strain to be a 'plane strain compression g in the X\- 
direction. Then, a numerical study of the variations of the quantity /3 + y/crf 
reveals that there exists a certain value v m such that for v > u m , there always 
exist values of A where this quantity is zero. This delimits a curve in the 
(A, z/)-plane beyond which ellipticity is lost. We find that this scenario is 
repeated for the equi-biaxial compression ^ in the x\- and ^-directions, 
and for the uni-axial compression Q in the Xi-direction. In particular, we 
find that v m ~ 0.272, 0.415, and 0.770, for the deformations @, Q, and 
(|4]), respectively, see Figure [2j 

3 Surface stability 

Now take a semi-infinite solid with strain-energy density ([T]). In this section, 
we study surface instability, and thus consider that the solid occupies the 
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x 2 > region in its deformed state, with the boundary x 2 = free of 
traction. 

Biot [1963] shows that when a half-space is subjected to a large homoge- 
neous deformation, there eventually exists a configuration where it is possible 
to superpose an infinity of small-amplitude static perturbations, with ampli- 
tude variations confined to the neighborhood of boundary. He considers that 
this configuration corresponds to the surface buckling, or surface instability, 
phenomenon. Dowaikh and Ogden [1990] find that, for a general isotropic 
incompressible solid, the surface stability criterion can be put in the form 

7>0, a- 7 + 2v^77 (/3 + 7 ) > 0. (8) 

Clearly, in view of the expression ^ 2 for 7, the first inequality is always 
satisfied for solids with strain energy ([I]), as is the case for neo-Hookean 
solids. Moreover, we recall that [Dowaikh and Ogden, 1990], if the second 
inequality is satisfied here, then the second inequality in ^ is automatically 
satisfied. In other words, surface instability always precedes internal (bulk) 
instability. 

We use criterion (j8j) 2 as an equality to compute the critical stretch of com- 
pression A cr for surface instability, in the cases of plane strain, as described 
by ([5]); of equi-biaxial compression, as described by ([3]); and of uni-axial com- 
pression, as described by Q. Figure [3] shows the resulting plots of A cr as a 
function of the LSS parameter v. At v = 0, we recover the critical stretch ra- 
tios of the neo-Hookean solid [Biot, 1963], i.e. 0.44 in uni-axial compression, 
0.54 in plane strain, and 0.66 in equi-biaxial compression. Then, there is a 
narrow range of values for ^ ~ where the half-space is slightly more stable 
than a neo-Hookean half-space. However, when v is greater than 0.15, we find 
that the half-space becomes unstable at higher compression ratios (earlier) 
than in the neo-Hookean case. This situation is in complete contrast with the 
isotropic models usually employed to model strain-stiffening biological soft 
tissues, which are always much more stable than the neo-Hookean solid with 
respect to surface buckling (see Goriely et al. [2006] and references therein). 

These results, and those of the previous section, suggest that tissues mod- 
eled with strain-energy density (jl]) become easily unstable when compressed. 
At the same time, it must be kept in mind that an infinite or semi-infinite 
medium does not exist and is an idealization; in Nature, all solids have finite 
dimensions. In the next sections, we investigate the influence of finite size 
for some prototype stability problems. 
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Figure 3: Surface instability of a compressed solid with strain-energy density 
given by (JTJ. For uni-axial compression, plane strain compression, and equi- 
biaxial compression, the surface instability occurs earlier than for the neo- 
Hookean solid (which corresponds to v — 0), once v > 0.15. 
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4 Slab stability 



Here we consider the homogeneous compression of a thick slab until it buckles 
(incrementally). Among others, Biot [1963b], Ogden and Roxburgh [1993], 
and Beatty and Pan [1998] have determined the bifurcation criterion in the 
case of a material with neo-Hookean strain-energy density. It reads 

tanh (A1A2 kh) 
tajoh(kh) 

when Ai is the principal stretch in the direction of compression. Here the plus 
(minus) sign is associated with antisymmetric/flexural (symmetric/extensional) 
modes of buckling, and kh gives a measure of the current slenderness of the 
slab. Explicitly, kh = mrh/l, where n is an integer, 2h is the current width of 
the slab (in the x\ direction), and 21 is its current height (in the x 2 direction). 
In terms of the slab's original width 2H and length 2L (say), we have 

kh = mr(\ 2 /\ 1 )(H/L). (10) 

Ogden and Roxburgh [1993] also give explicit bifurcation criteria for 
a general hyperelastic material, but there are several of these depending 
on whether the buckling is antisymmetric or symmetric and on the nature 
(purely imaginary or complex) of the roots to the following biquadratic: 

7 g 4 + 2/V + a = 0. (11) 

For neo-Hookean solids, these roots are qi i2 = ±i, g 3 4 = iiAiA^ 1 and their 
nature is fixed. For solids with strain-energy density Q however, the nature 
of the roots is strongly dependent of the value of the physical parameter v 
and of the geometrical parameters Ai, A 2 . Instead of discussing the different 
possibilities, we turn to the Stroh formalism to produce a numerically sta- 
ble bifurcation criterion, independent of the nature of the roots and of the 
buckling mode. 

First, we recall that the incremental equations of equilibrium can be put 
in the following form: 

ri{x 2 )' = ikNr)(x 2 ), (12) 

where rj is the displacement-traction vector and N the Stroh matrix. Ex- 
plicitly, taking the mechanical displacement u and the nominal traction s on 
the faces x 2 = const, in the form 

{u,s} = {U(x 2 ),ikS(x 2 )}e ikx \ (13) 



1 + Aj A 2 
4AA2 1 



(9) 
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where U and S are single-valued functions, we define 77 as if = [U, S]*. On 
the other hand, N is given, in the case where the x 2 = rkh faces of the slab 
are free of traction, by [Destrade et ai, 2005] 



N 



-1 I/7 

-1 

■(2/3 + 7) 0-1 

a-7 -1 



(14) 



Because N is a constant matrix, the general solution to (12) is of the 
form T](x 2 ) = N{x 2 )b, where b is a constant vector and M is a fundamental 
matrix solution: 



N(x 2 ) = [C\C 2 \C 3 \C 4 ] Diag(e lkqiX2 ,e [kq2X2 ,e [kg3X2 ,e ikq4X2 ) 



(15) 



with the g's and £' s eigenvalues and eigenvectors, respectively, of the eigen- 
problem N£ = qC,. Explicitly, the g's are the four roots of the characteristic 



equation det(iV— ql) = 0, which is the quartic (11 ), and the £'s are [Destrade, 
2007] 



C J = H?i/7, qj/l, -9i(g| - 1), q) - a 



(16) 



We may now write the solution on one face of the slab (at x 2 = h) in 
terms of the solution on the other face (at x 2 = —h) as: r](h) = Mr](—h), 
with M, the matricant [Shuvalov, 2000], defined as 



M = Af(h) [Af(-h)\ 



-1 



Mi 

M 3 M 4 



(17) 



where Mi, etc. are 2x2 submatrices. The incremental boundary conditions 
are that the end faces of the slab stay in sliding contact with the thrust 



platens: this is ensured by (10); and that its lateral faces remain free of 
traction: S(h) = S(—h) = 0. These identities lead to the robust form of the 
bifurcation criterion, 

and the left hand side of this equation is always real [Shuvalov, 2000] prior 
to buckling, independently of the nature of the roots q. 

We implemented this criterion for the case treated by Biot [1963] that is, 
plane strain compression ([2]), when the slab is not allowed to expand in the 
X3-direction. In Figure |4| we compare the critical stretches of compression for 
slabs made of solids with strain-energy density ((!]) (solid plots) with those for 
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Figure 4: Buckling of a compressed slab in plane strain: critical stretch 
against H/L, the slab's initial slenderness. Solid curves: solids with strain- 
energy density given by ([TJ. Dashed curves: neo-Hookean solids. 



neo-Hookean slabs (dashed plots). In the latter case we obtain two curves, 
the upper one corresponding to antisymmetric (flexural) modes of buckling 
and the lower one to a symmetric (compressional) mode, which is never 
attained, because it lies below the plots for antisymmetric buckling modes. 
In the former case, there is a single curve (solid line) for each value of the 
LSS parameter v, corresponding to the single bifurcation criterion (18). 

We plotted the value of the critical stretch of compression A cr against the 
initial slenderness of the slab H/L, for the range < H/L < 0.8, which is 
likely to cover all realistic dimensions of living tissues. We show the plots 
obtained at n = 1, but not those for n > 2, because they are all below the 
n — 1 antisymmetric mode, and are thus never reached. We see that in 
the infinitely thick slab limit H/L — > oo, the curves tend to the asymptotic 
value given by the critical stretch of surface stability of the previous section. 
In that limit, solid slabs with strain-energy given by Q buckle earlier than 
neo-Hookean slabs when v > 0.15, in agreement with the results of Figure [3} 
However soft tissues are in general quite slender and we note that in the low 
H/L region (0 < H/L < 0.2), the differences between the different curves 
become minute. In fact, the solids with strain-energy density ([I]) turn out to 
be more stable than neo-Hookean solids in the 0.2 < H/L < 0.3 range. 
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5 Tube stability 



Now we study the stability of a compressed tube. When it is placed between 
two lubricated platens, an elastic tube deforms homogeneously (in an equi- 
biaxial manner), and an exact incremental solution exists, due to Wilkes 
[1955]. 

The resolution of the stability problem can be conducted along lines sim- 
ilar to those outlined in the previous section, as shown by Goriely et al. 
[2008]. The main differences are that now the displacement and traction are 
sought in the form: 

{u,s} = {U(r),ikS(r)}e ike , (19) 

where r, 9 are the cylindrical coordinates in the deformed configuration; the 
displacement-traction vector is defined as rj = [U(r), rS(r)Y and it has six 
components instead of four; and the cylindrical counterpart to the rectangu- 



lar incremental equations of equilibrium (12) are 



%L(r) = -N(r)r,(r), (20) 
dr r 

where the 6x6 Stroh matrix now has variable components (see Goriely et 
al. [2008] for their explicit expression). The eigenvectors of N are written in 
terms of Bessel functions. The resulting bifurcation criterion is again of the 



form (18) 



The main tubular structures in the body are arteries, and we focus on 
this class of geometrical objects. We find in Delfino et al. [1997] that in their 
unloaded configuration, human carotid arteries have a wall thickness ranging 
from 0.4 to 0.9 mm, and inner diameter ranging from 1.8 to 3.3 mm. For our 
calculations we thus take B/A = 1.25 as being a representative geometric 
parameter, where A and B are the initial inner and outer radii, respectively. 
We then implement the numerical procedure of Goriely et al. [2008] to de- 
termine the critical stretch of compression for tubes with slenderness B/L 
in the 0.0 — 0.8 range, where L is the initial tube length. We find that in 
that range, the bifurcation curves corresponding to solids with LSS param- 
eter v = 0.25, 0.5, 0.9, all superpose almost exactly with the bifurcation 
curve of the neo-Hookean tube (at v = 0.0) and are indistinguishable one 
from another. For slender tubes, 0.0 < B/L < 0.21, the buckling occurs in 
the antisymmetric n = 1 mode. For stubbier tubes, 0.21 < B/L < 0.8, the 
barreling mode n = 2 is preferred. The curves of all other modes (n = 
and n > 3) are all situated below these two. Practically, it means that the 
tube cannot be compressed axially by more than 13% (because A cr = 0.87 at 
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Figure 5: Buckling of a compressed tube with dimensions close to those of 
a human carotid artery: B/A = 1.25 and 0.0 < B/L < 0.8, where A, B, L 
are the initial inner radius, outer radius, and length of the tube. Here we 
traced the bifurcation plots for v = 0.0 (neo-Hookean), 0.25, 0.5, 0.9, and 
found that all four superposed almost exactly. The figures on the left show 
the buckling modes corresponding to tubes with slenderness 0.1 (top figure) 
and 0.3 (bottom figure), denoted by two black spots on the dispersion curve. 

B/L = 0.21), at least when its slenderness B/L is less than 0.8. At that low 
level of strain, the solid with strain energy behaves almost exactly like a 
neo-Hookean solid, as can be checked by inspection of the curves in Figure 
[TJ This explains the coincidence of the curves in Figure [5j Only in the low 
part of the n = 1 curve can a distinction be made, but it is never expressed 
because the n = 2 buckling mode is prevalent at the corresponding ratios. 

6 Bending stability 

As a final example of compression instability, we study the effect of LSS 
on the bending stability criterion. For neo-Hookean solids, this instability 
occurs at compression ratios very close to the critical compression ratio of 
surface instability in plane strain. Quantitatively, the critical value for the 
circumferential stretch on the inner face of a bent neo-Hookean block is A cr ~ 
0.56, which corresponds to only 2% less compression than the surface stability 
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critical ratio of 0.54, see Section 3. At these levels of strain, the behavior 
of solids with strain-energy is likely to diverge from that of neo-Hookean 
solids. 

Take a rectangular block, made of a hyperelastic solid initially contained 
in the region 

< X 1 < L, -A<X 2 <A, 0<X 3 <H, (21) 
say. Then bend it into the cylindrical sector located in the region 

r a <r <r b , -a < 9 < a, < z < H, (22) 

where (r, 9, z) are the cylindrical coordinates in the current configuration; a, 
the angle of bending, is prescribed; and the inner and outer radii, r a and r&, 
respectively, can be found from Rivlin's bending deformation [Rivlin, 1949]: 



'2LX/a + (r2 + 7f)/2, 
Destrade et al. [2010] show that 



9 = aY/L, 



Z. 



2 2 




4A2 + ^T 2A 



(23) 



(24) 



Calling A2 the circumferential stretch, Destrade et al. [2010] also establish 
that it varies through the wall thickness from value Af at r = r a to value A2 
at r = r b , given by 

A^ = yV 1 + {^AafVf - (2Aa/L), X b 2 = + {^AafVf + (2Aa/L) 

(25) 

Now consider the possibility of a superposed incremental solution, de- 
scribing the apparition of prismatic wrinkles on the inner face of the bent 
block. Destrade et al. [2009] show that the corresponding incremental equa- 
tions of equilibrium can be put in the following fully non-dimensional Stroh 
form: 



-ry(A 2 



-Ar(A2)T7(A 2 ), 



(26) 



dA 2 A 2 

where the incremental displacement-traction vector 77 and the Stroh matrix 
iV are seen as functions of A2 only. Explicitly, 



N 



i -n 

— n(l — a/a) — i(l — a/a) —1/a 
«4i 1K12 — i — n(l — a/a) 

-i«i2 «22 —n i(l — a/a) 



(27) 
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where 

Ku = 2(/3 + a — cr) + n 2 [7 — (a — cr) 2 /a] , 
«i2 = n (2/3 + a + 7 - a 2 / a) , 

k 22 = 7- (a-cr) 2 /a + 2n 2 (/3 + a-cr), (28) 

and the quantities 7, a, and /3 are given in (J©J) , by taking the plane strain 
configuration (pi). Also, cr is the radial component of the Cauchy stress 
necessary to maintain the block in the bent configuration, with inner and 
outer faces free of traction: 

a = W(l/\ 2 , A 2 , 1) - W(l/\ a 2 , X a 2 , 1), (29) 

and n is the circumferential number, determined from the condition that 
there are no incremental normal tractions on the end faces 9 = ±a: it is 
such that [Haughton, 1999] 

n = pn/a, (30) 
for some integer p, the mode number. 



The differential system (26) must be integrated from A 2 to A 2 , subject to 
the incremental boundary conditions of traction-free bent faces: 

r 1 {\ a 2 ) = [U a M\ t,{)$ = [U b , Of, (31) 

where U a , Ub are non-zero constant vectors. The numerical treatment of this 
two-point boundary value problem is delicate due to potential stiffness issues. 
However, the compound matrix method (see [Haughton, 1999; Destrade et 
al. 2009]) smoothens these out. 

We implemented this technique to find the critical circumferential stretch 
ratio in terms of the block's aspect ratio L/A, see Figure [6J At 1/ = 0.0, 
the neo-Hookean solid buckles at A cr ~ 0.56, as previously shown [Haughton, 
1999], with a mode number which increases rapidly as the block is more 
slender (L/A increases). For v = 0.4, we find that blocks with strain energy 
([T]) buckle earlier, at A cr ~ 0.59. Similarly, A cr increases as v increases. Hence 
A cr ~ 0.61 at v — 0.8 and A cr ~ 0.62 at v — 1.2. Also, we see that there 
are less wrinkles on the inner bent face than in the neo-Hookean case. For 
example, a block which is 1.5 times longer than it is deep (L/(2A) = 1.5, 
giving L/A = 3) buckles with p = 4 wrinkles when v > 0.4, in contrast with 
p = 7 wrinkles in the neo-Hookean case [Destrade et al., 2009]. 

In conclusion, the level of compression at which instability occurs for the 
neo-Hookean solid is a good indicator of instability levels for solids with strain 
energy (§. Ifit occurs at high compression levels (as in this section), then the 
latter solids are unstable much earlier than neo-Hookean solids. Otherwise, 
(as in Sections [4] and |5| , there are no noticeable differences between the two 
types of solids. 
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Figure 6: Bending instability for a rectangular block with length-to-depth 
ratio L/(2A) G [0.25,1.75]. The neo-Hookean solid {y = 0.0) can be bent 
more than solids with strain energy density ([TJ, see plots at v = 0.4, 0.8, 1.2. 

Acknowledgments 

The first author gratefully acknowledges the support of a Senior Marie Curie 
Fellowship from the European Commission and the hospitality of Department 
of Continuum Mechanics and Structures at the E.T.S. Ing. Caminos, Canales 
y Puertos, Universidad Politecnica de Madrid. 

References 

[1] Beatty, M.F. [1996] "Introduction to Nonlinear Elasticity", In: Nonlin- 
ear Effects in Fluids and Solids, edited by M.M. Carroll and M. Hayes 
(Plenum Press, New York). 

[2] Beatty, M.F. and Pan, F. [1998] "Stability of an internally constrained, 
hyperelastic slab," International Journal of Non-Linear Mechanics 33, 
867-906. 

[3] Biot, M.A. [1963] "Surface instability of rubber in compression," Applied 
Science Research A 12, 168-182. 



16 



[4] Biot, M.A. [1963b] "Exact theory of buckling of a thick slab," Applied 
Science Research A 12, 183-198. 

[5] Biot, M.A. [1965] Mechanics of Incremental Deformations (John Wiley, 
New York). 

[6] Delfino, A., Stergiopoulos, N., Moore, J.E., and Meister, J.- J. [1997] 
"Residual strain effects on the stress field in a thick wall finite element 
model of the human carotid bifurcation," Journal of Biomechanics 30, 
777-786. 

[7] Destrade, M. [2007] "Interface waves in pre-stressed incompressible 
solids," In: Waves in Nonlinear Pre-Stressed Materials (M. Destrade 
and G. Saccomandi, Editors). CISM Lecture Notes (Springer, New York) 
495, 61-100. 

[8] Destrade, M., Gilchrist, M.D., Murphy, J.G. [2010] "Onset of non- 
linearity in the elastic bending of beams," ASME Journal of Applied 
Mechanics to appear. 

[9] Destrade, M., Ni Annaidh, A., and Coman, CD. [2009] "Bending in- 
stabilities of soft biological tissues," International Journal of Solids and 
Structures 46, 4322-4330. 

[10] Destrade, M., Ottenio, M., Pichugin, A.V., and Ro gerson, G.A. [2005] 
"Non-principal surface waves in deformed incompressible materials," In- 
ternational Journal of Engineering Science 42, 1092-1106. 

[11] Dowaikh, M.A. and Ogden, R.W. [1990] "On surface waves and defor- 
mations in a pre-stressed incompressible elastic solid," IMA Journal of 
Applied Mathematics 44, 261-284. 

[12] Gent, A.N. and Cho, I.S. [1999] "Surface instabilities in compressed or 
bent rubber blocks," Rubber Chemistry and Technology 72, 253-262. 

[13] Goriely, A., Destrade, M., and Ben Amar, M. [2006] "Instabilities in 
elastomers and soft tissues," Quarterly Journal of Mechanics and Ap- 
plied Mathematics 59, 615-630. 

[14] Goriely, A., Vandiver, R., Destrade, M. [2008] "Nonlinear Euler buck- 
ling," Proceedings of the Royal Society A 464, 3003-3019. 

[15] Han, H.-C. [2007] "A biomechanical model of artery buckling," Journal 
of Biomechanics 40, 3672-3678. 



17 



[16] Haughton, D.M. [1999] "Flexure and compression of incompressible elas- 
tic plates," International Journal of Engineering Science 37, 1693-1708. 

[17] Haughton, D.M. and Merodio, J. [2009] "The elasticity of arterial tissue 
affected by Marfan's syndrome," Mechanics Research Communications 
36, 659-668. 

[18] Eberth, J.F., Taucer, A.I., Wilson, E., and Humphrey, J.D. [2009] "Me- 
chanics of carotid arteries in a mouse model of Marfan syndrome," An- 
nals of Biomedical Engineering 37, 1093-1104. 

[19] Holzapfel, G.A., Gasser, T.C., and Ogden, R.W. [2000] "A new consti- 
tutive framework for arterial wall mechanics and a comparative study 
of material models," Journal of Elasticity 61, 1-48. 

[20] Merodio, J. and Haughton, D.M. [2010] "Bifurcation of thick-walled 
cylindrical sheLSS and the mechanical response of arterial tissue affected 
by Marfan's syndrome," Mechanics Research Communications 37, 1-6. 

[21] Minami, H. Asada, T. Gan, K., Abe, K. and Izumi, S. [2007] "My- 
ocardial ischemia due to compression of an unruptured thoracic aortic 
aneurysm in a patient with Marfan syndrome," General Thoracic Car- 
diovascular Surgery 2007, 248-251. 

[22] Mizuguchi, T. and Matsumoto, N. [2007] "Recent progress in genetics of 
Marfan syndrome and Marfan-associated disorders," Journal of Human 
Genetics 52, 1-12. 

[23] Ogden, R.W. and Roxburgh, D.G. [1993] "The effect of pre-stress on 
the vibration and stability of elastic plates," International Journal of 
Engineering Sciences 31, 1611-1639. 

[24] Rivlin, R.S. [1949] "Large elastic deformations of isotropic materials V: 
The problem of flexure," Proceedings of the Royal Society of London A 
195, 463-473. 

[25] Rodriguez, J. and Merodio, J. [submitted] "A new derivation of the 
bifurcation conditions of inflated circular cylinders of elastic material 
under axial loading. Application to aneurysm formation." 

[26] Shuvalov, A.L. [2000] "On the theory of wave propagation in anisotropic 
plates," Proceedings of the Royal Society: Mathematical, Physical and 
Engineering Sciences 456, 2197-2222. 



18 



[27] Watton, P.N., Heil, M., and Hill, N.A. [2004] "A mathematical model 
for the growth of the abdominal aortic aneurysm," Biomechanics and 
modeling in Mechanobiology 3, 98-113. 

[28] Watton, P.N. and Hill, N.A. [2009] "Evolving mechanical properties of 
a model of abdominal aortic aneurysm," Biomechanics and modeling in 
Mechanobiology 8, 25-42. 

[29] Wilkes, EW. [1955] "On the stability of a circular tube under end 
thrust," Quarterly Journal of Mechanics and applied Mathematics 8, 
88-100. 



19 



